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The two-flavour colour-superconducting (2SC) phase in small strangelets is studied. In order to 
describe the 2SC phase we use the three-flavour Nambu- Jona-Lasinio model. We explicitly take 
into account finite size effects by making use of the approximation for the density of states in 
spherical cavities called multiple reflection expansion (MRE) . The thermodynamic potential for the 
2SC strangelets is derived in the mean-field approximation with the help of the MRE. We found 
that 2SC phase survives in small strangelets with a sizable gap. Consequences for the 2SC phases 
are also discussed. 



I. INTRODUCTION 
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■ It is widely accepted that sufficiently cold and dense quark matter is a colour superconductor P, ■ At asymptoti- 
l/^ ' cally high densities the perturbative one-gluon exchange interaction, directly based on the first principle of QCD, has 
T-H , been used to clarify the properties of colour superconductivity. At moderate densities, on the other hand, the studies 

of QCD-motivated effective theories have revealed rich phase structures. 
, The most likely place for colour superconductivity is in the interior of compact stars. Therefore, the phases of quark 

■ matter with appropriate conditions for the interior of compact stars (i.e. neutrality conditions and /3-equilibrium 
conditions) has attracted a great deal of interest. The most striking feature of neutral and /3-equilibrated colour- 

■f—i • superconducting phases is the appearance of gapless colour superconductivity. In the case of a large strange quark 

■ mass, it has been revealed that the phase called gapless two-flavour colour-superconducting phase (g2SC) could 

■ appear as a ground state. Similar analyses have been done for the case of three flavours and clarified that the gapless 
colour-flavour- locked (CFL) phase could also appear in neutral and /3-equilibrated strange quark matter |^|^. (Note 
that, however, the gapless color superconductors develop a chromomagnetic instability 

f~i There exists other candidates for a colour superconductor called strangelets that are small chunks of strange quark 

Qh, matter ^ J^]. Strangelets with baryon number A <^ 10^ are free from (local) electric neutrality condition in contrast 
(JL) • to strange quark matter in bulk. As a consequence, the phases of strangelets could be different from that of bulk 
' strange quark matter. Instead, we need to take into account flnite-size effects explicitly in the study of strangelets. The 
properties of strangelets have been extensively studied by using various models with finite-size effects. In particular, 
^ ■ Madsen has shown that colour-flavour-locked strangelets are signiflcantly more stable than unpaired (normal) 
• ^ , strangelets, although the overall energy scale is still an open question because the bag constant is a phenomenological 
' input parameter. In previous preliminary work |l2l |. we investigated the behaviour of chiral symmetry in flnite size 
M , quark droplets consisting of up and down quarks within the two-flavour Nambu- Jona-Lasinio (NJL) model T^. In 
the NJL model, by contrast, the bag constant is generated d yna mically [l^ . To take into account flnite-size effects, we 
used the so-called multiple reflection expansion (MRE) [iM lla IitI Il8l| . In the MRE framework, the finite-size effects 
are included in terms of the density of states. The MRE has been used to calculate the thermodynamic quantities of 
finite-size quark droplets and reproduced well the results of mode-filling calculations ,18i] . 

In this paper, we focus on the quark core of strangelets with baryon number A < 10^ which are small enough to 
neglect electrons. In order to describe colour-superconductivity we use the three-flavour NJL model and study the 
conventional two-flavour colour superconducting (2SC) phase. Finite-size effects are incorporated in the thermody- 
namic potential by making use of the MRE. We deal with 2SC strangelets embedded in the physical vacuum. We 
do not discuss absolute stability of 2SC strangelets (as compared to a gas of ^^Fe), but how the 2SC gap behaves in 
strangelets and its effects on properties of strangelets. 

This paper is organised as follows. In Sec. II we formulate the thermodynamic potential for spherical 2SC 
strangelets. We concentrate on 2SC strangelets embedded in the physical vacuum and derive a set of coupled equa- 
tions so as to flnd stable (i.e. pressure balanced) 2SC strangelets. In Sec. Ill we restrict ourselves to zero temperature 
and present numerical results. Section IV is devoted to summary and discussions. 
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II. THERMODYNAMIC POTENTIAL FOR 2SC STRANGELETS 



To describe the 2SC phase we use the U(3)l x U(3)i?^ symmetric NJL model. We assmrie that the strange quark is 
sufhciently heavy and does not take part in Cooper pairing. The Lagrangian density is given by 

8 

£ - qiil'-d^ ~ m)q + ^ [{qr^qf + [qil^r'qf] 

1=0 

+Gd {qil5T^\^Cq^) {qCi-f5T^X\) (1) 

where q denotes a quark field with three fiavours {Nf = 3) and three colours {Nf. — 3), m — diag(0,0,TOs) is the 
quark mass matrix and C is the charge conjugation matrix, defined by C~^^^C = ~7j — ~C . The coupling 

constants Gs and Go have the dimension of (mass)"^. The Gell-Mann matrices t* (i = 1, • • • , 8) with t'^ = ^2/31 / 
act in the flavour space and is the antisymmetric generator of SU(3)c. 

We choose the model parameters (the ultraviolet cutoff and the coupling constant) as follows: Ayv = 0.6 GeV and 
Gs = 6.42 GeV~^. The diquark coupling constant is set to Gd = 3Gs/4, unless stated otherwise. The value of the 
strange quark mass shall be explicitly indicated in the following figures. 

In this paper we work in the mean-field approximation. In the mean-field approximation, we obtain the following 
Hamiltonian density: 

:^MFA = qin -y + m)q+ I- {qi-/5T'^X'^Cq'^ + h.c.) 



4Gr 
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where A = 2GD{qilb'r^X?C(f) is the 2SC gap parameter. Using HmfAj we can straightforwardly calculate the 
thermodynamic potential of bulk quark matter. 

In this work, however, we need to incorporate finite-size effects into the thermodynamic potential. To this end, we 
use the density of states derived from the MRE [isl H^. ItR IT^ . In the MRE framework, the density of states for a 
spherical system is written as k'^ pmke / ['2.1:'^) , where pmre = PMKE{k,m, R) is given by 

PMKE = l + ^fs -]+jr7^fc -]+■■■ ■ (3) 



kR-' \mj {kR) 

Here m denotes the (Dirac) mass of quarks and R is the radius of the sphere. The functions fsik/m) and fcik/m) 
represent the surface and curvature contributions to the fermionic density of states in the spherical cavity, respectively. 
The ellipsis corresponds to higher order terms in 1/R, which are neglected throughout. The functional forms of 
fsik/m) and fcik/m) are given by 

fs {—] = ^ ( 1 - -arctan — ) , (4a) 
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It should be noted that the functional form of fc for an arbitrary quark mass has not been derived within the MRE 
framework. The functional form of Eq. Ij4b|l is the ansatz by Madsen jlSj. Note also that Eqs. H4a|) and Ij4b|l have 
the following massless limits: 

lim /s(fc/TO) = , lim /c(fc/TO) = -l/(247r2). (5) 

m — *0 m — >0 

Figure 1 shows the MRE density of states. One can see that the finite-size effects reduce the density of states and 
Pmre becomes negative at small momenta. To avoid the unphysical negative density of states we shall introduce an 
infrared cutoff Air in momentum space. For the case of m 7^ 0, we numerically solve the equation pmre = with 
respect to k and use the larger one as Air. On the other hand, for the case of m = 0, the MRE density of states takes 
the form: 

PMRE(fc,-R) = 1 - 2(^kJi)2 ■ (6) 



FIG. 1: The density of states fc'^PMRE for m ^ 0. We numerically solve the equation pmre ~ and adopt the right-most root 
as infrared cutoff Ajr. 

Then, we obtain Air = \/2/{2R). We have confirmed that Air, behaves hke Air oc 1/R in both cases. 

Using the density of states we express the effective potential (per unit volume) of the spherical strangelets as 
follows 

fiMRE(A;Ai,r) = -^-2/ {e_ +e+ + 2rin(l + e-'^^-)(l + e-''^+)} 
4Ct_D Jmre 

-4 f |sgn(e_)i;_+£;+ + 2Tln(l + e-'3'^gn(^-)^-)(l + e-'5^+)| 
Jmre 

-3/ {£_+e+ + 2rin(l + e-'^^-)(l + e-''^+)}, (7) 

JMRE 



where e± = kzL^ are the quasiparticle energies of the unpaired (blue up and blue down) quarks, E± = ^J(k± /x)^ + 
are those of the gapped (red up, red down, green up, and green down) quarks, and e± — y^k"^ + ml ± /i are those of 
strange quarks. We have also introduced the sign function sgn(a;), sgn(a;) = ±1 for a; > or a: < 0, and the following 
short-hand notation. 



MRE 



-;r^PMRE- (8) 



Notice that we used a common chemical potential /i for all nine quarks. In general, the matrix of quark chemical 
potentials for three-flavour quark systems can be written as 

where Q, T3 and Tg are generators of U(l)em, U(l)3 and U(l)8, respectively. The indices i,j and a, b refer to flavour 
and colour, respectively. Hence, we should introduce three quark chemical potentials /ie, fJ-s and /xg in the study of 
strange quark matter in bulk. However, we can neglect the (local) electric charge neutrality and set fig ~ hi the 
study of small strangelets. We also neglect /i3 and /ig, for simplicity. 

Now we derive the gap equation with the MRE. For computations in a finite system, we temporarily choose a fixed 
radius R. Then, the extremum condition of f^MRE with respect to A reads 

A = 8Gd f [1 - 2Np{E^)] + ^ [1 - 2Nf{E+)]] . (10) 

Jmre l^- ^+ J 



This equation has a trivial solution (A = 0) as well as a nontrivial solution (A 7^ 0). The former corresponds to 
unpaired (normal) strangelets and the latter to colour-superconducting strangelets. 
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In order to look at stable (pressure balanced) strangelets, we employ the usual pressure balance relation between 
(the inside of) a strangelet and the outer physical vacuum. The inside pressure P can be defined by 



P = -Pmre - B, (11) 



where Pmre = -f^MRE and 



4Gs 



where M„ and Ms are the dynamically generated quark masses and they are the nontrivial solutions to the following 
equations, 

Of course, these equations can be derived from the extremum condition of the effective potential at T = ^ = 0. In 
Eq. ((TT|l . we have introduced the vacuum pressure (or bag constant) B so as to measure the pressure relative to the 
outer vacuum. We emphasise that we do not need to introduce the bag constant by hand. In the NJL model, the bag 
constant is generated dynamically. 

We first solve the gap equation [Eq. (|10|) ] and pressure balance relation P ~ [Eq. (|ll|l ] self-consistently. Then, 
the baryon number A of the stable strangelet is computed by using the following relation 

A = VnB, (14) 

where V = AttR^ /3 is the volume of the spherical strangelet and ns is the baryon number density. The baryon number 
density, which is one third of the quark number density Uq, is obtained by taking the partial derivative of JImre: 

ns = "-g/3, (15) 

d^MRE 

= 4 / (Npie-) - NFie+) - [1 - 2iVj.(i?_)] + |^ [1 - 2Nf{E+)] + ^ [iVj.(£_) - Np{e+)]] . (16) 

JMRE I Z ) 

where Nf{x) = l/{e^^ + 1) is the Fermi distribution function. 

Having solved the above three equations self-consistently, we can find stable colour-superconducting strangelets (i.e. 
A dependence of A, /i and R). We can also compute other thermodynamic quantities. The most important, perhaps, 
is the energy per baryon number E/A of strangelets: 



E 
A 



_ _£_ 



(17) 



where the energy density £ is given by 



^MRE + ^ fJ-fnf + B. (18) 

f—u,d^s 



Here nf is the number density of quark flavour / that satisfy the relation X]f ~ "-g- virtue of the common 
chemical potential n and the pressure balance relation, one can see that the following relation holds: 

E , , 

- = 3m. (19) 

^ p=o 
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III. NUMERICAL RESULTS 



In this section we focus on strangelets at zero temperature and present numerical results. Before proceeding to the 
results, we give a brief survey of the finite-size effects. Suppose a finite-size quark droplet embedded in the physical 
vacuum. It is generally known that the finite-size quark droplet is energetically disfavoured as compared to quark 
matter in bulk. In other words, the finite-size effects increase the quark chemical potential. Hence, our results show 
a size dependence (i.e. baryon number dependence). 

Let us start with the baryon number dependence of the 2SC gap (see Fig. 2). We present the results for the cases 
771s = and TTis = 0.12 GeV, for comparison. It is clear that, for relatively large baryon numbers {A > 100), the 
gap remains approximately constant. The result indicates that finite-size effects becomes less important for large 
baryon numbers (i.e. large radii), as they should. The difference in the size of the gap between the case rus — and 
771s = 0.12 GeV is the result of the difference in the chemical potentials of the gapped quarks. At fixed A, the effect 
of rris is to increase the chemical potentials of up and down quarks. Consequently, the gap grows with 77is. 

In contrast, the curves behave differently at small baryon numbers {A < 100). In the case of ttIs = 0.12 GeV, as A 
is decreased the gap decreases. The decrease of the gap is caused by the fact that the chemical potentials of gapped 
quarks are increased by finite-size effects and, then, they are close to the UV cutoff (see Fig. 3). 

In the case of ttIs — 0, the behaviour of the gap at small baryon numbers is also due to the finite-size effects. 
However, the chemical potentials of the gapped quarks are not close to the UV cutoff (Fig. 3). As A is decreased, the 
chemical potentials of gapped quarks are increased by the finite-size effects. Hence the gap grows with decreasing A. 
As A is further decreased the gap reaches a peak and, finally, drops. The decrease of the gap at very small baryon 
numbers {A < 10) arises from the fact that the density of states near the Fermi surface is decreased by the IR cutoff. 

Figure 4 shows energy per baryon number E/A of 2SC strangelets as a function of A. The strange quark masses are 
set to 771s — 0.12 GeV. We also present the results for the cases of strong diquark couphng {Go = Gs) and unpaired 
strangelets. It appears that 2SC strangelets are more stable than unpaired strangelets and the strong diquark coupling 
significantly stabilises 2SC strangelets. The results can be easily understood as follows. In the 2SC phase, the pairing 
energy contribution to the thermodynamic potential is roughly given by (volume term only) 

nfc) ^ _ AV ^ (20) 

that is the product of the pairing energy and the number density of gapped quarks near the Fermi surface [T^ . 
Therefore, it is quite reasonable that 2SC strangelets are more stable than unpaired strangelets. It is also natural 
that the strong diquark coupling (i.e. large paring energy) makes 2SC strangelets energetically more favoured. For 
small A the energy per baryon number increases dramatically because of the growth of the chemical potential by the 
finite-size effects. 

To examine E/A of various phases we compute E/A for unpaired, 2SC, and CFL strangelets as a function of A. 
(See Fig. 5. The curve for CFL strangelets is taken from Ref. "2^). The diquark coupling and the strange quark 
mass are taken to be G_d = 3Gs/4 and ttIs = 0, respectively. The pairing energy in the CFL phase takes the following 
form [T^ : 

Hence the pairing energy contribution to the thermodynamic potential in the CFL phase is larger than that in the 
2SC phase because all nine quarks participate in Cooper pairing. As one would expect, CFL strangelets are more 
stable than 2SC strangelets. In general, however, we need to include nonvanishing ?7is and study the competition of 
2SC and CFL phases in detail. 

It is interesting to look at the strangeness fraction of 2SC strangelets in the light of their production in laboratories. 
The strangeness fraction is given by 

/. = "\ ■ (22) 

nu + 71<J + 71s 

The fraction for unpaired strangelets should satisfy /s = 1/3 at ?7is = 0. Figure 6 shows the strangeness fraction of 
2SC strangelets for the cases of ttIs = 0, ttIs = 0.12 GeV, and ttIs = 0.15 GeV. The diquark coupling is taken to be 
Gd = 3G5/4. The fraction deviates from /s = 1/3 even in the case of ms — 0. This deviation can be understood as 
follows. The number density of an unpaired quark is given by 

no = 2 / 9ifi- k), (23) 
Jmre 
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FIG. 2: The 2SC gap A as a function of baryon number A for the cases of tUs — and nis = 0.12 GeV. The diquark couphng 
is taken to be Gd = 3Gs/4. 

where 0{fJ-— k) is the occupation number for a noninteracting massless quark. On the other hand, the number density 
of a gapped quark is given by 



where + /i — fc)/(2_E_) is the occupation number of a gapped quark. One can easily check that the relation 
nA > no holds. Hence, the number of gapped quarks increases in the 2SC phase and the relation fs < 1/3 holds 
even in the case of TOs = 0. For the cases of massive strange quark, the fraction decreases further because the Fermi 
momentum of strange quarks, kp — ^/i^ — m^, decreases due to the mass gap. One may argue that fs is reduced 
further if one takes into account the dynamically generated mass of the strange quark. However, it has been shown 
that the finite-size effects enhance the restoration of chiral symmetry fl^. Therefore, at present, the effect of the 
dynamical strange quark mass is an open question. The strangeness fraction in CFL strangelets is, presumably, larger 
than that of 2SC strangelets, because CFL pairing enforces the equal Fermi momenta for all quarks. The complete 
calculation including dynamical quark masses is an interesting problem to be investigated in future work. 

Figure 7 shows the densities achieved in 2SC strangelets as a function of A. 2SC strangelets with massive (m^ = 0.12 
GeV) strange quarks have higher densities than that with massless strange quarks. This is understood by considering 
that nonzero TOj gives a contribution to the pressure. The leading-order contribution from nonvanishing rus is roughly 
given by (volume term only) 



Comparing the cases of TOs = and iris — 0.12 GeV, we have confirmed that fls overpowers the small difference of 

(2SC) 

flp . Then 2SC strangelets with massive strange quarks require larger chemical potentials to maintain the pressure 
balance and, therefore, have larger density than that with massless strange quarks. 

Inversely, the radii of strangelets behave as shown in Fig. 8. The size of 2SC strangelets at very small baryon 
numbers {A ~ 10) is around 1-2 fm. The size is comparable with the coherence length of the Cooper pair ^ = l/(7rA). 
Then, it is uncertain whether the Cooper pair exists in such a small system or not. In these regions, mode-filling 
calculations may be available ■ 




(24) 



IV. SUMMARY AND DISCUSSIONS 



In summary, we studied the properties of two-flavour colour-superconducting strangelets. We used the Nambu- 
Jona-Lasinio model to describe colour superconductivity and the multiple reflection expansion to take account of 
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FIG. 3: The quark chemical potentials /i of the pressure balanced 2SC strangelets as a function of baryon number A. Recall 
that the UV cutoff is taken to be Ayv = 0-6 GeV. 
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FIG. 4: The energy per baryon number E/A as a function of A for 2SC strangelets (solid lines) and unpaired strangelets 
(dotted line). 

finite-size effects. We formulated the tfiermodynamic potential for 2SC strangelets including finite-size effects in 
terms of the MRE density of states. We then solved the set of coupled equations; the gap equation, the pressure 
balance relation and the baryon number condition. Wc found that a sizable 2SC gap survives in small strangelets 
[A ~ 100) though its behaviour at small baryon numbers depends on the strange quark mass. We also found that 
2SC strangelets are more stable than unpaired strangelets due to the pairing energy. It should be noted here that 
2SC strangelets are not absolutely stable as compared to a gas of ^^Fe. We did not examine the competition with 
other possible phases (e.g. colour-flavour-locked phase, chirally broken phase, and so on). We would like to make a 
comment on the sensitivity of the results to model parameters. Using several sets of parameters, we examined the 
sensitivity of our results. We found that the results do not change much with the choice of Ayv and G5, as long 
as they are fixed by fitting physical quantities at T = /x = 0. Rather than Ayv and G5, the diquark coupling has 
an effect on the results. A large diquark coupling yields a large 2SC gap and reduces the strangeness fraction. For 
instance, the strangeness fraction decreases to ~ 0.25 at = 0.12 GeV and Gd/Gs = 1-2 due to the large 2SC 
gap (A ~ 0.28 GeV). However, such a strong coupling may induce the fluctuation of pairing field and, then, invalidate 
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FIG. 5: The energy per baryon number EjA as a function of A for CFL strangelets (bottom), 2SC strangelets (middle) and 
unpaired strangelets (top). The diquark coupling and the quark mass are set to Gd = 3Gs/4 and nis = 0, respectively. 
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FIG. 6: The strangeness fraction fs of 2SC strangelets as a function of baryon number A. The diquark coupling is set to 
Gd = 3Gs/4:. The strange quark masses are taken to be rus = 0, 0.12, 0.15 GeV (from top to bottom). 

the present description of the 2SC phase. 

The MRE has problems concerning its reliabihty. First, Pmre should have terms proportional to l/R^ and 

so on. They arc dominant at small radii. Further, the MRE causes a negative density of states. Although we simply 
avoided the latter problem by introducing an infrared cutoff, these problems should be solved in the future. However, 
for relatively large systems {A > 100), these problems have rather minor effects and the results presented in this 
paper would hold. 

Finally, we comment on the outlook for future studies. It is very interesting to study the competition with other 
phases, taking account of the realistic value of the strange quark mass. In particular, inclusion of the chirally broken 
phase and the colour-flavour-locked phase would affect the present analysis. It would be also interesting to examine 
colour-superconducting strangelets embedded in the vacuum at finite temperature and/or density. 
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FIG. 7: The baryon number density tib in 2SC strangelets divided by the normal density no = 0.17fm 
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FIG. 8: The radius R of 2SC strangelets as a function of baryon number A for rus = 0,0.12 GeV. Both curves behave like 
R oc A^/^ 



Acknowledgments 

I would like to thank D. Rischke and I. Shovkovy for discussions and a critical reading of the manuscript. 



[1] B. Barrois, Nucl. Phys. B129, 390 (1977); S.C. Frautschi, in it Proceedings of the Workshop on Hadronic Matter at 

Extreme Energy Density, Ericc, Italy, 1978, edited by N. Cabibbo and L. Sertorio (Plenum, New York 1980); D. Bailin 

and A. Love, Phys. Rop. 107, 325 (1984). 
[2] K. Rajagopal and F. Wilczek, in At the Frontier of Particle Physics - Handbook of QCD: Boris loffe Festschrift, edited 

by M. Shifman (World Scientific, Singapore, 2001); M. Alford, Annu. Rev. Nucl. Part. Sci. 51, 131 (2001); D.H. Rischke, 

Prog. Part. Nucl. Phys. 52, 197 (2004). 
[3] I. Shovkovy and M. Huang, Phys. Lett. B 564, 205 (2003). 
[4] M. Huang and I. Shovkovy, Nucl. Phys. A729, 835 (2003). 



10 



[5] M. Alford, C. Kouvaris, and K. Rajagopal, Phys. Rev. Lett. 92, 222001 (2004). 
[6] M. Alford, C. Kouvaris, and K. Rajagopal, Phys. Rev. D 71, 054009 (2005). 
[7] M. Huang and I. Shovkovy, Phys. Rev. D 70, 051501 (R) (2004); 70, 094030 (2004). 
[8] R. Casalbuoni, R. Gatto, M. Mannarelh, G. NarduUi, and M. Raggieri, Phys. Lett. B 605, 362 (2005). 
[9] See, for example, J. Madsen, astro-ph/9809032 C. Greiner and J. Schaffner-Bielich, nucl-th/9801062i 
[10] E. Farhi and R.L. Jaffe, Phys. Rev. D 30, 2379 (1984); M.S. Berger and R.L. Jaffe, Phys. Rev. C 35, 213 (1987); 44, 
566(E) (1991); C. Greiner, D.H. Rischke, Horst Stocker, and P. Koch, Phys. Rev. D 38, 2797 (1988); G. Greiner and H. 
Stocker, Phys. Rev. D 44, 3517 (1991); J. Madsen, Phys. Rev. Lett. 70, 391 (1993); E.P. Gilson and R.L. Jaffe, Phys. Rev. 
Lett. 71, 332 (1993); J. Madsen, Phys. Rev. D 47, 5156 (1993). 
[11] J. Madsen, Phys. Rev. Lett. 87, 172003 (2001). 
[12] O. Kiriyama and A. Hosaka, Phys. Rev. D 67, 085010 (2003). 
[13] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961); 124, 246 (1961). 

[14] K. Schertler, S. Leupold, J. Schaffner-Bielich, Phys. Rev. C 60, 025801 (1999); M. Buballa and M. Oertel, Phys. Lett. B 

457, 261 (1999). 
[15] R. Balian and C. Bloch, Ann. Phys. (N.Y.) 60, 401 (1970). 
[16] E. Farhi and R.L. Jaffe, Phys. Rev. D 30, 2379 (1984). 

[17] M.S. Berger and R.L. Jaffe, Phys. Rev. G 35, 213 (1987); 44, 566(E) (1991). 
[18] J. Madsen, Phys. Rev. D 50, 3328 (1994). 

[19] M. Alford and K. Rajagopal, J. High Energy Phys. 06, 031 (2002). 
[20] O. Kiriyama, hep-ph/0401075 

[21] P. Amore, M.G. Birse, J. A. McGovern, and N.R. Walet, Phys. Rev. D 65, 074005 (2002) 



